Critical exponents of the disorder-driven superfiuid-insulator transition in 
one-dimensional Bose-Einstein condensates 
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We investigate the nature of the superfiuid-insulator quantum phase transition driven by disorder 
for non-interacting ultracold atoms on one-dimensional lattices. We consider two different cases: 
Anderson-type disorder, with local energies randomly distributed, and pseudo-disorder due to a 
potential incommensurate with the lattice, which is usually called the Aubry- Andre model. A scaling 
analysis of numerical data for the superfiuid fraction for different lattice sizes allows us to determine 
quantum critical exponents characterizing the disorder-driven superfiuid-insulator transition. We 
also briefly discuss the effect of interactions close to the non-interacting quantum critical point of 
the Aubry-Andre model. 
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I. INTRODUCTION 



A superfiuid-insulator transition in a disordered nonin- 
teracting system of bosons at zero temperature is a spe- 
cial type of quantum phase transition (QPT) [![. Instead 
of the more conventional competition between different 
interactions, it is disorder that causes a drastic change in 
the nature of the ground state, thus altering the physical 
characteristics of the material. A similar type of tran- 
sition from a metal to an insulator, usually called the 
Anderson localization transition, was first proposed by 
Anderson 0, and has been extensively studied in elec- 
tronic systems In general, the approach focus on 
the conductance behavior as the Fermi level changes in 
the vicinity of the mobility edge, which separates local- 
ized and extended one-particle states. In this context, 
the lower critical dimension has been determined to be 
dr, = 2, which means that all the states are localized in 
one dimension for any finite amount of disorder. Never- 
theless, given that the states in a strictly non disordered 
system are extended, there is a clear change of regime 
when the disordered strength is reduced to zero, which 
can be characterized as a QPT. 

In the past decade, enormous progress in the tech- 
niques for creating ultracold atom systems in laboratory 
settings extended the interest in the disorder effects and 
localization to bosonic systems (for recent reviews, see 
[3-0] ). For bosons, the transition is from the insulator 
(localized) state to the superfiuid one. It was observed 
both for laser-speckle disorder Q and quasiperiodic opti- 
cal lattices [9( in Bose-Einstein condensates of 87 Rb and 
39 K atoms, respectively. While speckle disorder comes 
close to the Anderson-type random disorder considered 
in theoretical approaches, quasi-periodic lattices present 
a superposition of the lattice potential with an incom- 
mensurate one and can be viewed as experimental real- 
izations of pseudo-disorder models like the Aubry-Andre 
(AA) model [lo| . The latter also shows superfiuid and 
localized regimes in one dimension, but the transition be- 



tween them occurs at a nonzero critical disorder (lCj, lll( . 

Recently, we have investigated numerically the 
superfiuid-insulator transition in one-dimensional, nonin- 
teracting systems of bosons with these two types of disor- 
der [l2| . Here we focus on the scaling properties of the su- 
perfiuid fraction near the superfiuid-insulator transition, 
obtaining the relevant critical exponents. For random 
disorder, even though the superfiuid phase is destroyed 
for arbitrarily weak disorder, we show that the transition 
can still be described as a quantum critical phenomenon 
with well-defined critical exponents and power-law scal- 
ing behavior. The same happens for the A A model, but 
the universality classes arc different. 

Our starting point is a well-known scaling relation for 
the singular part of the superfiuid density p s close to a 
quantum superfiuid-insulator phase transition (l3| . 
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where g measures the distance to the quantum critical 
point (QCP), v is the correlation length exponent (i.e., 
the correlation length diverges as £ ~ at the QCP), 
d is the spatial dimension, and z is the dynamic critical 
exponent associated with the QCP. The superfiuid den- 
sity is directly related to the helicity modulus [Til ], and 
can be viewed as a measure of the system's response to a 
phase-twisting field. Thus, it is natural to interpret the 
correlation length as a phase-coherence length. In the 
insulating phase it should coincide with the localization 
length, which measures the spatial extent of the wave 
functions. This holds also for disordered metals Q. 

In a finite system, even at criticality, the correlation 
length is limited by the system size L, and the finite- 
size-scaling form of the superfiuid density is 
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The corresponding relation for the superfiuid fraction 
(f s = L d p s ) is 



f s ~L-^F(L\gn- 



(3) 
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This last equation is suitable to determine the critical 
exponents v and z from a numerical evaluation of f s for 
various lattice sizes, as we do in the following. 



II. ANDERSON-LIKE DISORDER 

The usual Hamiltonian describing interacting bosons 
on a lattice is known as the Bose-Hubbard Hamiltonian, 
and is given by 



H 



]eini+Qy^(alaj+a\ai)+—y^ni(ni-l) , (4) 



where a\ and a i are the creation and annihilation oper- 
ators of a boson at the lattice site i, rii — a\a i is the 
corresponding number operator, each site has a single 
bound state of energy Si, hopping between sites is re- 
stricted to nearest neighbors, with amplitude O, and U 
is a local repulsive interaction. In the rest of this paper, 
energies are measured in units of the tunneling amplitude 
f2. An Anderson-like disorder [2j is introduced by choos- 
ing random local energies with a uniform distribution in 
the range — A/2 < £j < A/2, so that A is a measure of 
the disorder strength. 

We carried out a thorough numerical analysis of the 
above model for the non-interacting case in one spatial 
dimension. Details of this numerical study are given in 
Ref. [12[. We recall one of the main results reported 
there, namely that the superfiuid fraction for a lattice of 
size L obeys the relation f s = cxp(— A/A^) 4 / 3 , where 
Al is a characteristic disorder strength for suppression 
of superfluidity, which scales with the lattice size as 
Al = CI -3 ' 2 . This latter relation is consistent with the 
expected value of the critical disorder strength A c = 
for destroying the superfiuid phase in the thermodynamic 
limit for a one-dimensional system. Furthermore, defin- 
ing g = A — A c = A, this scaling of A^ with L is rec- 
ognized as the finite-size version of the general relation 
£ ~ Isl - "' immediately yielding the correlation-length 
exponent v — 2/3. 

Equation Q implies that L z f s is a universal func- 
tion of LA". The corresponding plot of our data for 
different lattice sizes is shown in Fig. [TJ where it is clear 
that all the data collapse onto a universal curve. The 
appearance of f s alone as the scaling quantity in the 
vertical axis means that the dynamic critical exponent 
is z — 2. The scaled variable of the horizontal axis in 
the collapsed plot confirms the value v = 2/3 for the 
correlation-length exponent. Actually, for the present 
problem, we are able to determine explicitly the scal- 
ing function F(x) in Eq. ([3]). The above mentioned ex- 
pression for f s implies that F(x) = cxp(— x 2 /C 2 ). Since 
x = L/£ = LA U , the value z = 2 for the dynamic ex- 
ponent implies a jump of the superfiuid fraction from 
to 1 at the transition occurring for A — in the limit 
L — > 00. This jump is reminiscent of that of the helicity 
modulus in the two-dimensional XY model 11511 . 
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FIG. 1: (Color online) Finite size scaling of the superfiuid 
fraction for Anderson- type disorder. According to Eq. Q, the 
horizontal-axis variable for collapsed curves is L\g\ v , which 
gives v = 2/3, while the absence of any rescaling of f s implies 
that z — 2. 



The value v = 2/3 obtained here for the correlation- 
length exponent has not been determined previously, to 
the best of our knowledge. This new exponent for the 
supcrfluid-insulator transition seems to violate the in- 
equality v > 2/d, which holds for other disordered sys- 
tems [l3|. However, this inequality has been proved only 
for interacting systems and for nonzero critical values of 
the parameter driving the transition, which is not the 
case here. On the other hand, the dynamic exponent 
z = 2 implies that the effective dimension of the quan- 
tum phase transition [l| is d e s = d + z = 3. For dis- 
ordered interacting bosons, the Bose-glass-to-superfluid 
transition is characterized by the relation z = d [13l |. 
Therefore, interacting and noninteracting bosons are in 
different universality classes with respect to the localiza- 
tion transition. In the renormalization group language, 
interaction is a relevant term close to the disordered non- 
interacting fixed point. 



III. AUBRY-ANDRE MODEL 

The Aubry- Andre model can also be described by the 
Hamiltonian with U — 0, except that the distribution 
of local energies is not random, but periodic with a period 
incommensurate with the lattice spacing. These energies 
are usually written as 



Acos(27r/3i), 



(5) 



where (3 = (1 + \/h)/2 is the golden ratio and i assumes 
integer values from 1 to £. This is actually a special 
case of the Harper model [16] for electrons in a two- 
dimensional lattice in the presence of a perpendicular 
magnetic field, for which Eq. ([5]) holds for any value of 
/3, with different characteristics of the spectrum for ra- 
tional or irrational values. Disorder-like effects here are a 
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FIG. 2: (Color online) Finite size scaling of the superfluid 
fraction for the AA model. The data collapse in different 
curves for even and odd numbers of lattice sites (respectively, 
lower and upper curves). 
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consequence of the incommensurability between the "ex- 
ternal potential" and the lattice. Aubry and Andre [Toj 
proved that for this model localization occurs only when 
the strength of the potential A is larger than the criti- 
cal value A c = 2. For finite lattices, it is convenient to 
replace (3 with f3 n = F n+ i/F n , the ratio of two consec- 
utive Fibonacci numbers, whose limit for n — > oo is the 
golden ratio. Then, the lattice size must be chosen as 
L = F n in order to allow for the use of periodic bound- 
ary conditions. For this kind of finite lattices, the critical 
value A c = 2 remains a rigorous result [Tl| , since it corre- 
sponds to a duality between the Hamiltonians in position 
and momentum space. 

It was shown in Ref. [l2[ that the superfluid fraction 
undergoes a very sharp transition around A = 2 for es- 
sentially all lattice sizes. This sharpness makes it difficult 
to directly extract the correlation-length exponent, as 
done for random disorder (Sec. [TTJ). Here, we concentrate 
on a narrow region around A c , searching for the appropri- 
ate scaling variable proportional to g = A — A c = A — 2, 
and the appropriate scaling of the superfluid fraction. 
Our results are shown in Fig. [3J The data collapse onto 
two universal curves, for even and odd numbers of lat- 
tice sites. Although the scaling functions are different 
for these two cases, the critical exponents for which the 
curves collapse are the same. In view of Eq. ([3]) we im- 
mediately identify the correlation-length exponent v = 1 
from the x-axis scaling variable in Fig. [2] and the dy- 
namic exponent z = 2.374 from the y-axis scaling. It 
was already known [n| that v = 1 for this model. Next 
we discuss the obtained value of z in the light of proper- 
ties of the energy spectrum. 

The spectrum of the Harper model has been thor- 
oughly studied in the past [l7H2l| . For general rational 
values of (3 it is multifractal at A = 2, yielding the fa- 
mous Hofstadter butterfly [22j , shown in Fig. [3] There 
we highlight the case that we are studying here, with 
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FIG. 3: (Color online) Spectra of the Harper model, high- 
lighting the spectrum corresponding to the AA model (for 
a rational approximation of the golden ratio /3 — 987/610). 
Its fractal nature is illustrated in the bottom by expanding 
the small box drawn inside the middle panel. We show the 
spectrum for /? — 1, which is the same as for /3 according to 
Eq. ©. 



/3 being a rational approximant of the golden ratio. In 
particular, the figure shows our result obtained from nu- 
merical diagonalization of the Hamiltonian for a lattice 
of size L = 610. The two bottom panels illustrate the 
fractal nature of this spectrum. 

With our replacement of /3 by a ratio of two Fibonacci 
numbers, /3 n — F n+ \/F n , the spectrum is equivalent 
to the one for /3„ = J3 n — 1 = F„_i/F„, which con- 
tains F n bands and F„_i gaps. As discussed in de- 
tail in Refs. [ill when F n = L — > oo the width 
AEl of a given band belonging to the spectrum scales 
as AEl ~ L~ 7 , with different regions of the spectrum 
associated with different values of 7 (not to be confused 
with the susceptibility critical exponent). In particular, 
a maximum value 7 max = 2.374 corresponds to band- 
edge states. On the other hand, the band width is a 
characteristic energy of the system and therefore should 
scale as £ _z , which means that AEl ~ L~ z . Our find- 
ing of z = 7 max is in agreement with the relevant state 
for the zero-temperature superfluid-insulator transition 
being the bottom edge of the lowest-lying band. 
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FIG. 4: (Color online) Zero-temperature phase diagram near 
the localization transition in the presence of a small interac- 
tion. The color scale indicates values of the superfluid frac- 
tion. The line U c (g) in the vicinity of the non-interacting 
QCP is essentially linear (with a slope close to 0.1) implying 
that it is dominated by the analytic part f(g) in Eq. (J5J). 



IV. EFFECTS OF INTERACTION 

The interaction term in Eq. ((4} can be treated as a 
relevant field close to the QCP. The knowledge of the 
dynamic exponent allows us to generalize the scaling re- 
lations close to the QCP for small but finite U. The free 
energy, for example, is given by 



F s ^\g\^ d+ ^P(U/\ 9 r). 



(6) 



where again we used the fact that U is an energy, and thus 
scales with £~ 2 . From the above equation, we see that 
the scaling contribution to the critical line separating the 
superfluid and insulating phases is 



U c (g) oc \g\ 



(J) 



The critical exponents are those associated with the dis- 
ordered non-interacting QCP at U = 0, A = Ac = 2. 
Since the product vz — 2.374 is large we have to take 
into account analytic contributions to the shape of the 
critical line. This line can be written in general as 



U c = f(g) + a±\gl 



(8) 



where f(g) is an analytic function, and ± refers to the 
sign of g. Expanding close to the QCP, analytic contri- 
butions up to the second order dominate over the scaling 
term when g — > 0. To illustrate this point, in Fig. [4] we 
show a phase diagram close to the noninteracting fixed 
point at A c = 2, for repulsive (U > 0) and attractive 
(U < 0) interactions. This is a plot of the superfluid 
fraction (color scale) as a function of A and U, obtained 



by diagonalizing the Hamiltonian ((4]) for N — 8 interact- 
ing bosons on a lattice of L = 8 sites. We can see that a 
straight line [i.e., f(g) ~ g] is a very good approximation 
to the boundary between the superfluid and localized re- 
gions. Even though the transition is smoothed out by the 
small lattice size, it is worth mentioning that the value 
of f s at the critical point (U = 0, A = 2) is compatible 
with the lower curve of Fig. [2] for L = 8. 



V. CONCLUSIONS 

We studied the superfluid-insulator transition for 
bosons on a one-dimensional lattice, both with random 
disorder and the pseudo disorder described by the Aubry- 
Andre model, the two prototype models employed in the 
investigation of localization for ultracold atoms in optical 
lattices. Using a finite-size-scaling analysis of the super- 
fluid fraction, we obtained the critical exponents char- 
acterizing this transition. The superfluid fraction yields 
the correlation-length exponent v and the dynamic crit- 
ical exponent z. For random disorder we found v = 2/3 
and z — 2, while for the AA model the results are 
v = 1 and z = 2.374. The other critical exponents can 
be obtained from the quantum hyperscaling relations [l[ 
2 - a = f(d + z) and 2/3 = v{d + z - 2 + rf). These 
two models fall into different universality classes, which 
is not surprising since the critical disorder strength for 
the superfluid-insulator transition is zero for Anderson- 
like disorder and nonzero for the AA model, which also 
exhibits a multifractal energy spectrum at the QCP. 

It is interesting to observe that the scaling form of 
the free energy for nonzero temperature T can be used 
to determine the thermodynamic behavior close to the 
superfluid-insulating QCP. From it, a general dependence 
of the specific heat with temperature is obtained [l3| . 
with the form C ~ T d ' z . 

Our brief discussion of interaction effects in the AA 
model shows that the critical point moves to stronger 
disorder for repulsive interaction, and to weaker disor- 
der in the attractive case. Regions of the phase diagram 
that correspond to localized and superfluid regimes in 
the thermodynamic limit are separated by a line that is 
approximately linear, reflecting the dominance of nonsin- 
gular contributions. 
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